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The observation of isolated positive and negative charges, but not isolated magnetic north and 
south poles, is an old puzzle. Instead, evidence of effective magnetic monopoles has been found 
in the abstract momentum space. Apart from Hall-related effects, few observable consequences of 
these abstract monopoles are known. Here, we show that it is possible to manipulate the monopoles 
by external magnetic fields and probe them by universal conductance fluctuation (UCF) measure- 
ments in ferromagnets with strong spin-orbit coupling. The observed fluctuations are not noise, but 
reproducible quasiperiodic oscillations as a function of magnetisation direction, a novel Berry phase 
fingerprint of the magnetic monopoles. 



Quantum states in solids are classified by a crystal mo- 
mentum vector and a band index. The space spanned by 
the momentum vectors is known as the momentum space. 
Each band index defines an energy-band of allowed elec- 
tronic energy levels in the momentum space. Momentum- 
space magnetic monopoles arise from energy-band cross- 
ings pQ. Each band crossing point produces a magnetic 
monopole with a quantised topological magnetic charge, 
characterised by a Chern number [1 . An electric parti- 
cle traversing a closed curve in momentum space accumu- 
lates a geometric phase from the monopole fields [H[3J . So 
far, these abstract monopoles have revealed themselves 
only through Hall-related effects [H [5] , but we show that 
they can also be manipulated and probed by UCF mea- 
surements. 

UCFs are observed experimentally as reproducible 
fluctuations in the conductance in response to an applied 
external magnetic field B The fluctuation pattern 
is known as the magneto-fingerprint of the sample [6]. 
Recent experiments on the ferromagnetic semiconductor 
(Ga, Mn)As report two different B c periods in the con- 
ductance fluctuations pj [5] : a slow, conventional oscilla- 
tion for high magnetic fields and a much faster oscilla- 
tion for low fields when the magnetisation rotates. The 
present work reinterprets these recent experimental re- 
sults and shows that the fast oscillations are caused by 
the relocation of momentum-space magnetic monopoles. 
Rotation of the magnetisation relocates the monopoles, 
which leads to a geometric phase change of the closed 
momentum-space curves. The numerical results demon- 
strate, in good agreement with the experiments, that a 
geometric phase change is observed with fast UCF oscil- 
lations, implying a novel Berry phase fingerprint of the 
monopoles. 

The underlying physics of UCFs is quantum interfer- 
ence between different paths across the sample [6] . Let A c 
denote the quantum mechanical probability amplitude 
for propagating along the classical path x c (£). The am- 
plitude can be expressed as A c = |A C | exp(iS[x c (t)]/h) in 
terms of the action S[x(t)] = JdtL(x, x), where L(x,x) 
is the Lagrangian, x = dx/dt, and |A C | is the prob- 



ability to follow the path x c (i) [5]. When an exter- 
nal magnetic field B is applied, the term — e J dtx • A 
should be added to the action [5], where e is minus 
the electron charge, and A is the vector potential cor- 
responding to B = V x A. Let us separate out the 
magnetic field-dependent phase and rewrite the ampli- 
tude as A c = A c exp(— ie/h J dtx c ■ A). The conductance 
G is proportional to the total probability of propagat- 
ing across the sample, G(B) oc E c ^c| 2 - Reformulat- 
ing the line integral associated with the vector potential 
as a surface integral using Green's theorem, one finds 
G(B) cx Ecc' A*A C , exp(i27r<f cc ,(B)/<I> ), where ^(B) 
is the magnetic flux enclosed by the loop formed by the 
paths x c (i) and x c ' (t) and <f> = h/e. Changing the mag- 
netic field randomises the phase difference between differ- 
ent pairs of paths, causing the conductance to fluctuate. 
A typical period B c of these quasiperiodic oscillations 
is when the dominant paths experience a relative phase 
shift of 2?:. Assuming that typical paths approximately 
enclose the sample area A, leads to B c = $ /^4 [5]. 

A closed loop in real space also corresponds to a closed 
loop in momentum space. In systems with either bro- 
ken inversion or time-reversal symmetry, there is also a 
phase associated with paths in momentum space [lj [3J . 
Semiclassically, this Berry phase effect is included in the 
Lagrangian as hA^ ■ k, where A'")(k) = i (u n |Vk| u n ) 
is the Berry connection, \u n ) is the periodic part of the 
Bloch function, and n is the band index [3J. The prop- 
agation amplitudes accumulate a geometric phase factor 
exp(i J dk ■ A^(k)) along a path in momentum space. 
A closed momentum-space curve acquires a phase equal 
to the flux of the effective field (k) = V k x A( n ) (k) 
that the loop encloses [U [3J . The effective field is known 
as the Berry curvature [U [2] : 

^ n (E n (k) - E m (k)) 2 

where H is the Hamiltonian of the system and E n (k) 
is the dispersion relation of the nth band. Momentum- 
space magnetic monopoles are singularities in the Berry 
curvature where energy bands cross at isolated points [4 , 
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5 . In ferromagnets with strong spin-orbit coupling, the 
Hamiltonian is not invariant under rotation of the mag- 
netisation [10 . Changing the magnetisation direction 
relocates the magnetic monopoles, inducing a geomet- 
ric phase change in the propagation amplitudes. The 
external magnetic field can rotate the magnetisation in 
UCF experiments on ferromagnets. The phase change 
of a closed real-space curve then also acquires impor- 
tant contributions from the geometric phase change of 
the corresponding closed momentum-space curve. We 
demonstrate that the magnetic monopoles give rise to 
fast conductance oscillations at low magnetic fields. This 
novel and large magnetic monopole effect is qualitatively 
different from the studies of Berry phase effects in two- 
dimensional electron gases with Rashba spin-orbit cou- 
pling since these systems exhibit no effective momentum- 
space monopoles [TT]. Also, the effect we compute quan- 
titively differs from the weak peak splitting effects seen 
therein by 1 order of magnitude. 

In the following discussion, the Berry phase effect on 
UCFs will be investigated for the ferromagnetic semi- 
conductor (Ga, Mn)As. The system is modeled by the 
Hamiltonian [TU] 
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(P- J) 2 + h- J + V(r). (2) 



The band structure of the host compound is described by 
the two first terms in Eq. ([2]), characterised by the Lut- 
tinger parameters 71 and 72 . J is a vector of 4x4 spin ma- 
trices for 3/2 spins, and P = p — eA is the canonical mo- 
mentum operator in the presence of an external magnetic 
field B = V x A. m e denotes the electron mass. The 
third term describes the exchange interaction between 
the holes and the local magnetic moments, modeled by 
a homogenous exchange field h. To model disorder, we 
used the impurity potential V(r) = J2i ViS(r— Rj), where 
Vt and Ri are the strength and position of impurity num- 
ber i and 5(r) is the delta function. 

The magneto crystalline anisotropy in (Ga, Mn)As is 
complicated and depends on several material parame- 
ters such as doping, strain and shape: see Ref. [TU] 
and references therein. We consider two cases: 1) a per- 
pendicular easy magnetisation axis that is valid, for ex- 
ample, for (Ga, Mn)As grown on (Ga, In)As and 2) a 
uniaxial in-plane easy magnetisation axis that is valid 
for (Ga, Mn)As bars grown on a GaAs substrate [TU] , 
Here, the magnetisation and hence h are assumed to 
be governed by the following magnetic free energy e = 
K u sin 2 (^/) - MBcos(4> M - <f> B ), where (f> M (4>b) is the 
angle between the exchange field h (applied magnetic 
field B) and the current. 

For the numerical UCF calculation, wc considered a 
discrete rectangular conductor sandwiched between two 
clean reservoirs with rectangular cross sections defined 
by L x = 30 nm and L z — 14 nm. The spacing between 
the lattice points is a x — a y = a z = 1 nm, significantly 



smaller than the typical Fermi wavelengths used at ~ 
5 nm. We assumed one impurity at each lattice site in 
the conductor. The current direction is [010], and the 
crystal growth direction and applied magnetic field are 
along [001]. We used the Landau gauge A = Bxy. For 
direct comparisons with experimental findings, we used 
parameters appropriate for (Ga, Mn)As: 71 = 7.0, 72 = 
2.5, Fermi energy E F = 78 meV and |h| = 31 meV. The 
impurity strengths Vi are uniformly distributed between 
— Vo/2 and Vq/2. Vq = 0.75 eV, which leads to a mean 
free path of ~ 6 nm. We assume that M = 2 x 10 4 A/m 
and the uniaxial anisotropy constant K u = 5 x 10 3 J/m 3 , 
giving the anisotropy field B u = 2K U /M — 0.5 T, which 
is similar to the experimental value found in Ref. [7]. 
The Landauer-Biitikker formula is used to calculate the 
conductance from a stable transfer matrix method |12j . 
More details about the numerical calculation method can 
be found in Ref. IT51. 
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FIG. 1: (a) Dispersion curves along the [001] axis when the 
exchange field h is pointing along [001]. Each band crossing 
point ±ki,2 gives rise to a momentum-space monopole. The 
black dotted line is the Fermi level used in the numerical UCF 
simulation, (b) Energy surfaces in the k x k y plane near the 
ki band crossing point in figure (a), (c) Energy surfaces in 
the k x k y plane near the k2 band crossing point in figure (a). 

Let us first analyse and classify the monopoles and 
then use simple semiclassical considerations to estimate 
the Berry phase-induced UCF oscillation period. The ge- 
ometric phase-induced conductance oscillations appear 
for weak external magnetic fields, and we can there- 
fore neglect the real-space magnetic field in the anal- 
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ysis of the Berry curvature. Without magnetic fields 
and disorder, the Hamiltonian in Eq. ^ has four band 

crossing points located at ±|ki| = ±y(|h| m e )/(27 2 ?i 2 ) 
and ±|k 2 | = ±\ 



|h| m e )/(72?i 2 ) along the momentum- 
space axis parallel to the exchange field, as shown in 
Fig. [Tk. Each crossing point gives rise to a magnetic 
monopole. The eigenfunctions of the Hamiltonian are of 
the form exp(ik • r) x n k where x n k i s a four-component 
spinor. Because the helicity operator £ = k • J/ |k| com- 
mutes with the Hamiltonian when k is parallel to h, 
the eigenspinors along this axis in momentum space are 
X mz k||h = U(9, 4>) \m z ) where \m z ) are eigenvectors of J z 
and U (0, (f>) is the unitary rotation operator that rotates 
the quantisation axis parallel to h. At the point ki , x 3 ; 2 
and X1/2 are degenerate eigenspinors. Close to this point, 
the Hamiltonian couples these two states only weakly to 
the X-3/2 an d X-i/2 states, and it can therefore be writ- 
ten as a 2 x 2 matrix in the basis of (1 0) T = X1/2 an d 
(0 1) T = X3/2- Expanding the Hamiltonian around the 
degenerate point, H = H (ki) +5k- Vk 1 H (5k = k — ki), 
treating the last term as a perturbation and consider- 
ing the case when h — |h| z, we obtain the local 2x2 
Hamiltonian: 



H = E (<5k) I + lx • <r, 



(3) 



where Eo (<5k) is an energy shift of the two energy 
bands away from the source point, er is a vector of 
Pauli matrices, I is the identity matrix, and x = 
272ft 2 |ki| /m e (—\/35k x , y/ZSky, 26k z ). When the crys- 
tal momentum k varies in time, the effective Hamiltonian 
in Eq. ([3| describes a spin in a time-varying magnetic 
field and the electron accumulates a well-known geomet- 
ric phase from the Berry curvature field [2] 



n (±) (k) = T 



2|k|' 



(4) 



where =F refer to the upper and lower energy bands near 
ki. We have here reparameterised the momentum space 
as k H> k = (— x\, x 2 , X3) for clarity. The topological 
magnetic charge of this monopole, its Chern number, is 

l/27r/ s2 ^ (±) -fidS-Tl EQ- 

Because the Berry curvature is inversely proportional 
to (E n (k) — E m (k)) 2 , bands that nearly cross over a 
larger region in momentum space produce a stronger 
monopole. The structures of the energy bands in 
Figs. [Tp-c show that the k 2 monopole is stronger than 
the ki monopole. A similar simple perturbative analysis 
of the k 2 monopole cannot be carried out, but numeri- 
cally, we find that the curvature decays asymptotically 
as k 2 , and have a Chern number of ±2. 

As can be seen from Fig. [Ik, the Berry curvature field 
from the two ±ki monopoles is experienced by orbits in 
the third and fourth bands, whereas the ±k 2 monopoles 



give a geometric phase effect only to momentum-space 
curves in the second and third bands. In the third band, 
the ±ki and ±k 2 monopoles have topological charges of 
opposite signs and therefore counteract each other. Or- 
bits in the lowest band do not experience any monopole 
field. Therefore, paths on the Fermi surface of the second 
band, which are experiencing the strong ±k 2 monopoles, 
dominate the Berry phase-induced conductance fluctua- 
tions. 
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FIG. 2: (a) The second-band Fermi surface and the Berry 
curvature field (black arrows) on this surface when the ex- 
change field points along [001]. A typical path is shown in (a) 
as a yellow curve where fc z ~ ir/L z . The curve accumulates a 
geometric phase of ~ 2n. (b) The second-band Fermi surface 
and the Berry curvature field (black arrows) on this surface 
when the exchange field is rotated by an angle of 0.13. The 
yellow curve represents the same real-space curve as in (a). It 
accumulates a vanishing geometric phase. In both plots, the 
red dot is the k2 monopole source, (c) The geometric phase 
change of the curve as a function of the rotation angle of the 
exchange field. The geometric phase angle and rotation angle 
are in units of, respectively, 2n and 7r/20 rad. 



A typical closed momentum-space curve on the second- 
band Fermi surface is shown in Fig. [2k. The correspond- 
ing real-space curve is found from the semiclassical equa- 
tion Tly = VkE 2 (k)|E, = E F . Because the Fermi surface is 
not rotationally symmetric, the momentum-space curve 
corresponding to this fixed real-space curve changes lo- 
cation on the Fermi surface when h is rotated, as illus- 
trated in Fig.[2]D. The associated geometric phase change 
of the relocated momentum-space curve, calculated nu- 
merically, is shown in Fig. [2];. We found that a rotation 
angle on the order of 0.13 changes the Berry phase of 
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this fixed real-space curve by 2it. The origin of the rapid 
phase change occurs when the curve in Fig. [2^, encloses 
the strong Berry curvature field region on the Fermi sur- 
face, shown as black arrows in Fig. [2^-b, whereas the 
curve in Fig. [2]d is relocated outside this region. 

The magnetic field needed to rotate the magnetisation 
by 0.13 rad is therefore an estimate of the oscillation 
period of the Berry phase fingerprint. Using the free en- 
ergy defined above, this leads to the oscillation period 
£jBerry ^ B u /10, where B u is the minimal external mag- 
netic field needed to align the magnetisation along the 
hard magnetic anisotropy axis. 
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FIG. 3: The conductance versus the applied magnetic field 
for the perpendicular magnetisation easy axis (black) and in- 
plane easy axis (red) for two different impurity configurations. 
The upper curve is shifted 1.5e 2 /h upwards for clarity. The 
grey region marks the region where the sample's magnetisa- 
tion changes from in-plane to perpendicular. The curves for 
the perpendicular and in-plane easy axes are coincident for 
B > B u = 0.5T because the same impurity configuration is 
used. Inset: a close-up for small magnetic fields. 

Let us next confirm our semiclassical analysis by a nu- 
merical UCF simulation of the system including disorder 
and magnetic field in Eq. (|2j). We consider two cases: 1) 
a perpendicular easy magnetisation axis that is valid, for 
example, for (Ga, Mn)As grown on (Ga, In) As and 2) a 
uniaxial in-plane easy magnetisation axis that is valid for 
(Ga, Mn)As bars grown on a GaAs substrate [10] . The 
external magnetic field is applied along the [001] growth 
direction. 

First, consider the case of a perpendicular magnetisa- 
tion easy axis where the magnetisation is aligned along 
the growth direction. We see in Fig. [3] that the conduc- 
tance has a weak increasing trend for increasing B. Here, 
the magnetic field squeezes the spatial extension of the 
wave function [14] . allowing more conducting channels 
to be open for increasing B. Imposed on the increasing 
trend, there are strong conductance fluctuations with a 
dominant period B c ~ IT. Because the magnetisation 
here is always along [001], the only change in the quan- 
tal phases comes from the magnetic flux. For a wire of 



A = 30 x lOOnm 2 xy area, the dominant fluctuation pe- 
riod is B c = $o/^4 ~ IT [6], consistent with the data 
shown in Fig. [3] 

Second, consider the case of an in-plane easy axis 
where h rotates from the [010] direction to the hard 
[001] axis when the magnetic field increases from to 
the anisotropy field B u = 0.5T. The decreasing trend 
of the conductance for increasing B 6 [0, B u ], shown in 
Fig. [3] is the standard anisotropic magnetoresistance 
effect [T2| . Imposed on the decreasing trend, the con- 
ductance fluctuates wildly for B < B u with a period on 
the order of i?„/10. Here, changing B leads to changes 
in the direction of h, which relocates the position of the 
momentum- space magnetic monopoles and thereby the 
geometric phase for a given real-space orbit. This gives 
rise to the extraordinarily fast conductance fluctuations 
shown in Fig. [3] For B > B u , the Berry phase is fixed 
and the UCF again exclusively comes from the conven- 
tional magnetic flux. 

Similar to what is found for the intrinsic anomalous 
Hall effect [H [5] , the effect is strongest for Fermi energies 
near the monopole sources. We expect the effect to also 
be present for more highly doped (Ga, Mn)As systems 
that require a six- or eight-band model in which more 
monopoles are expected to exist. 

In conclusion, the UCF simulation in Fig. [3] semiquan- 
titatively reproduces the experiments in Refs. 0[S], and 
together with our semiclassical analysis, it reinterprets 
the fast oscillations as a Berry phase fingerprint. 
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